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^ Abstract 

JL . We compute the spectrum of thermal photons created in Au+Au collisions at ^/snn = 200 

CD 

GeV, taking into account dissipative corrections in production processes corresponding to the 
quark-gluon plasma and hadronic phases. To describe the evolution of the fireball we use a viscous 
fluid dynamic model with different parametrizations for the temperature-dependence of w/s. We 
find that the spectrum significantly depends on the values of w/s in the QGP phase, and is almost 
insensitive to the values in the hadronic phase. We also compare the influence of the temperature- 
dependence of w/s on the spectrum of thermal photons to that of using different equations of state 
in the fluid dynamic simulations, finding that both effects are of the same order of magnitude. 
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I. INTRODUCTION 



There is currently great interest and effort in understanding the transport properties 
of the quark-gluon plasma (QGP) and of the hadronic matter created in ultrarelativistic 
heavy ion collisions at RHIC and LHC jl-3|. As evidenced by the large elliptic flow that is 
measured, matter created in these events behaves as a nearly perfect fluid with a very low 
viscosity-to-entropy ratio rj/s. For this reason, relativistic viscous fluid dynamics with rj/s < 
0.4 has been very successful in describing the hadronic observables such as the radial and 
elliptic flows [l[ 0^12 ] . For recent reviews on the application of relativistic hydrodynamics 
to heavy ion collisions see e.g. 

aa. 

Great efforts are currently focused on developing new theoretical tools to compute more 
accurately the tra nsp ort coefficients of the QGP and hadronic matter from microscopic 



models (see e.g. 
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18 n as well as on extracting them more precisely from RHIC and 



LHC measurements 2Ml2l|. Although it is expected that the transport coefficients of the 
QGP depend on temperature, the impact of a temperature-dependent i]/s on momentum 
anisotropics as obtained from simulations has been investigated only recently 
These studies show that the momentum anisotropics are indeed affected by the tempereature 
dependence of rj/s, pointing to the conclusion that further investigation is necessary to be 
able to extract precise values of rj/s from data. 

In this paper we focus on the spectrum of thermal photons produced during the expansion 
of the fireball created in Au+Au collisions at ^/snn = 200 GeV. In constrast to hadrons, 
the photons created in the interior of the fireball pass through it without any interaction, 
;hus giving information on the properties of bulk nuclear matter rather than on its surface 



24j-|37|- 



The main purpose of this work is to quantify the sensitivity of the spectrum of thermal 
photons to the temperature-dependence of rj/s. To this end, we compute the spectrum of 
thermal photons as obtained from fluid dynamic simulations using simple parametrizations 
of rj/s as a function of temperature, which allow us to disentangle the impact that the 
temperature dependence of rj/s in the different phases (QGP and hadronic) has on thermal 
photon spectra. 

Thermal photon spectra have been calculated within the framework of ideal hydrody- 



namics (see |25|), using the Israel-Stewart formalism 27H29| and more recently using a 
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viscous (3 + 1)D imp 
phase was studied in 



ementation 



27 



351 ] . The inclusion of viscous corrections during the QGP 



28j using Israel-Stewart theory and taking into account both shear 



and bulk viscosity, while in 29[ the nonequilibrium correction to photon production due 
to Compton scattering and qq annihilation were computed at leading-log order. Recently, 
thermal photon spectra including dissipative corrections has been computed in a (1 + 1)D 
hydrodynamic model including shear and bulk viscosities and with an EOS corresponding to 



a first order phase transition 



37j . The spectrum of thermal photons has also been computed 



using hydro- kinetic hybrid approaches to describe the evolution of the fireball 32M34| . In 



26l | we compared the thermal photon spectra obtained from second order fluid dynamics and 



a divergence-type theory that includes dissipative effects to all orders in velocity gradients 



12] , concluding that differences in dissipative fluid dynamic models are a significant source 



of uncertainty in the precise determination of rj/s from data. It is important to notice that 
there has been an increasing interest in the elliptic flow of photons produced in heavy ion 
collisions, since this observable may provide deeper understanding on the properties of the 
strongly coupled QGP - see [35|, |36[ where the elliptic flow and the spectrum of photons are 
calculated. 

It is appropriate at this point to comment on the limitated scope of our work. On the 
fluid dynamics side, we neglect bulk viscosity, which peaks at the critical temperature and is 



known to increase photon production 
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Besides, our model is (2 + l)D boost invariant. 
This approximation is not so limiting if one is concerned only with the dynamics near the 
mid-rapidity region, which is the case studied here. We also assume that the baryonic chem- 
ical potential is zero in the central rapidity region, and therefore the conservation equation 
of the net baryon number is not considered. Moreover, we use smooth initial conditions for 
the energy density. The importance of initial-state fluctuations and its impact on the hydro- 
dynamic flow and final observables has only recently been noted (see 38M40l| ) . In particular, 



as shown in 



41 



421 ]. event-by-event fluctuations of the initial energy density produced in 



heavy ion collisions at RHIC significantly enhance the production of thermal photons as 
compared to a smooth initial-state averaged profile in ideal hydrodynamic simulations. On 
the photon production side, we do not take into account prompt and jet-medium photons, 
whose contribution to the total photon spectrum becomes important at large transverse 



momentum 
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3l| . In spite of these limitations, the model employed here describes realis- 



tically the hydrodynamic evolution of matter created at RHIC (see [8|, [l2j, |23J]), allowing us 
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to quantify the impact of the temperature-dependence of rj/s on the spectrum of thermal 
photons. 

This paper is organized as follows. In Section [TT] we describe the viscous fluid dynamic 
model and provide numerical details of the simulations, and also describe the procedure 
used to compute the photon spectra. In Section (Till we present and discuss our results, and 
finally we conclude in Section IIVI 



II. THEORETICAL MODEL 



In this section, we first give a brief overview of the equations of viscous fluid dynamics 
used to model the evolution of matter created in heavy ion collisions, and provide details 
of the simulations performed. Then we go over to describe the processes that taken into 
account to compute the thermal photon spectra. 



A. Viscous fluid dynamics 

We will now briefly review the (2 + 1)D hydrodynamic model used in this paper. A 



detailed description of the model as applied to heavy ion collisions can be found in |8|, 112 ] . 

In what follows, Latin indices stand for transverse coordinates (x,y), is the geomet- 
ric covariant derivative, D = u^D^ and V M = A^ V D V are the comoving time and space 
derivatives, respectively, and brackets around indices imply taking the spatial, symmetric 
and traceless projection of a tensor. We employ Milne coordinates defined by proper time 



t = \Jt 2 — z 2 and rapidity ip = arctanh(z/t), and work in flat space-time. We assume boost 
invariance so that all quantities are independent of ip- The fluid velocity is u = (u T , u x , u y , 0) 
and is normalized as u^u^ = 1. The stress-energy tensor reads T^ u = pu^u v — pA^ v + IP", 
with A^ u = — u^u v ', where p and p are the energy density and the pressure in the local 
rest frame, and IP y is the viscous shear tensor. For a conformal fluid as considered here, 
Tjf = 0, so p = 3p and the bulk viscosity vanishes. 

The hydrodynamic equations are the conservation equations for the stress-energy tensor 
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together with the evolution equation for the shear tensor IT" 7 . The latter reads [16l . [1 



i a 



d T U ia = —U ia VX —IT* + — — — cr 

3u T Tt.U 1 ' T 7r U T 

1 : U< l U a> ^ ^ ^Du^ (1) 



2t 7T t] 2 u t m u t 

— a-rr 

where r\ is the shear viscosity, (r^Ai) are second-order transport coefficients, and = 
y<^ M ^> j g fix-st-order shear tensor. 

The conservation equations for T^ v and the evolution equation for IP" given in Eq. ([1]) 
are solved selfconsistently using the parameters defined in the following section. 



B. Numerical setup 

We choose (p, u x , u y , U xx , U xy , U yy ) as independent variables. Solution of the hydrody- 
namic equations requires initial conditions for the six independent variables, which we take 
to be u x = u y = and (IF 2 , U xy , U yy ) = 0, while the initial energy density profile is calcu- 



lated using a simple Glauber model 12]. As in previous fluid dynamic simulations of heavy 
ion collisions at ^snn = 200 GeV performed with the same model [3| , we take the initial- 
ization time to be To = 1 fm/c, use a computational grid of 13x13 fm and use values for the 
second-order transport coefficients corresponding to a strongly-coupled J\f = 4 Super- Yang 



Mills (SYM) plasma Q, Q ll6Hl8|] , ?v = 2(2 - ki2)rj/(sT) and Ai = 77/(2ttT), where s is the 
entropy density and T is the temperature. 

To solve the hydrodynamic equations it is necessary to provide the value of r;/s as an 
input, as well as an equation of state (EOS) relating the pressure and energy density. Here 
we will employ three parametrizations of 77/s as a function of temperature. Specifically, we 
consider the cases in which the ratio stays constant throughout the entire fluid dynamic 
evolution (Model A), or varies in the hadronic phase and stays constant in the QGP phase 
(Model B), or else varies in the QGP phase and stays constant in the hadronic phase (Model 
C). The three cases are shown in the upper panel of Figure [TJ We note that for Model A 
we choose the value of 77/s to be equal to the average value obtained from Models B and C, 
namely rj/s = 0.11, which gives sense to the comparison between the three parametrizations. 
This value of i]/s is on the lower side in terms of fitting charged hadron elliptic flow to data 



With respect to the EOS, we employ the EOS obtained by Laine and Schroder [431 ]. which 
connects a high-order weak-coupling perturbative QCD calculation at high temperatures to a 
hadron resonance gas at low tem per atures, via an analytic crossover as suggested by Lattice 



QCD calculations - see e.g. 



44j-|49| . We notice that this EOS is the same as that used in 



previous simulations focusing on charged hadron elliptic flow 



12] . The EOS is shown 



in the lower panel of Figure [TJ together with a different EOS 23J that is used later on for 
comparison (see Section ITTT|) . 

Since the model used here is purely hydrodynamic, a cut-off for the value of 77/s must be 
imposed in order to avoid the breakdown of the fluid description. The value of rj/s at which 
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FIG. 1: (Color online) Three different parametrizations for the temperature dependence of rj/s: 
constant (Model A), constant in the QGP phase and varying in the hadronic phase (Model B), and 
constant in the hadronic phase and varying in the QGP phase (Model C) (upper panel). Equation 
of state used in the simulations (denoted by LQCD) that was obtained by Laine and Schroder 43 , 
together with a different EOS (denoted by EOS 1) that is used for comparison (lower panel). 
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it is sensible to impose this cut-off is constrained by the comparison of results obtained from 
hydrodynamic and kinetic simulations to data. In our simulations we se t rjl s < 0.4, which 
is close to the upper bound imposed by such comparisons to data - see [23|. We note that 
the conclusions extracted from our results do not depend on the precise value of the cut-off 
for the value rj/s. 



C. Thermal photon production 



In order to compute the spectrum of thermal photons created during the evolution of the 
fireball we consider the processes of Compton scattering, qq annihilation and bremsstrahlung 
in the QGP phase, and tttt — > p7, np — > tt^ and p — > 7T7T7 in the hadron phase 



25 



26(. In all 



cases, the nonequilibrium distribution function of the quarks, /(a; M ,p M ), is calculated from 



Grad's quadratic ansatz 



12|: 



/(^) = / + /o(l-/o) 



(2) 



2T2(p + p) 

where fo is the equilibrium distribution function of the quarks. The evolution of the tem- 
perature T, energy density p, fluid velocity w M and shear tensor are obtained from the 
fluid dynamic model described in the previous section. 

The production rate for Compton scattering and qq annihilation reads 



25 



29] 



E 



dN 



-aa. 



^eJ)T 2 /(^,^)ln 



cE 
a.JT 



(3) 



with c ~ 0.23, a = 1/137 and [50j a s {T) = 6tt/((33 - 2JV» ln(8T/T c )), where N f is the 
number of quark flavors, e/ is the electric charge of the quark and T c is the critical temper- 
ature, which we set to be T c = 170 MeV in agreement with Lattice QCD simulations 44]. 
The contribution due to qq annihilation when an additional scattering is included is 

dN 8 



E- 



d 4 xd 3 p 37T 5 
with Jt ~ 1.11 and Jl ~ 1.06. 



aa, 



(J2e 2 f)ETf(x»,p»)[J T - J L ] 
f 



(4) 
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For bremsstrahlung we have 5l| 
dN 



E 



d 4 xd 3 p 



7T 



x 



/ 

3C(3) + ^ + ^ In 2 + 4Li 3 (-/(x^)) 
2 |Li 2 (-/(^,^)) - |^ln(l + 



(5) 



where ( is the zeta function and Li m = Y^nLi zU / nTn are P°lyl°g functions. 

For the production rate of photons durin g th e hadronic phase we use the estimate cal- 
culated by Steffen and Thoma |52| (see also 53|]), which reproduces the sum of production 
rates for the processes tttt — > p^, up — > and p — > titc-/: 

dN 



E- 



4 .8T 2 - 15 e- 1 /( 1 - 35£T ) f{x^pf 



(6) 



dx 



(7) 



d 4 xd?p 

The total thermal photon spectrum is then obtained by integrating the sum of the pro- 
duction rates over the evolution of the fireball as obtained from fluid dynamics, thus 

/ dN 
\d?p T dY 

where Tip are the initial and final times of each phase, Y is the rapidity of the nuclei, 
dx = (dx, dy) and pt is the transverse momentum. Note that the photon energy in the 
comoving frame is given by p^cosh(y — Y'). In what follows we will limit ourselves to the 
case Y = Y' = 0. 



- r* E dN 
Y d 4 xd 3 p 



III. RESULTS 

In this section we present and analyze the spectra of thermal photons computed from the 
fluid dynamic model described above. 

Figure |2] shows the contributions to the thermal photon spectrum corresponding to the 
QGP (upper panel) and the hadronic phases (lower panel), for Models A (temperature- 
independent), B (varying in the hadronic phase) and C (varying in the QGP phase) for rj/s. 
It is seen that the difference in the spectra corresponding to Models A and B is very small 
in the whole range of values of pt that we consider. In contrast, the spectrum computed 
from Model C is significantly larger, and the difference with the spectra of Models A and B 
increases with increasing values of pr- It is also seen that the difference between the results 
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obtained from Model C and those obtained from Models A and B is slightly larger for the 
hadronic phase. 

Previous studies focusing on the elliptic flow of charged hadrons computed from viscous 
fluid dynamics with different profiles of rj/s as functions of temperature showed 
that at RHIC energies the elliptic flow is almost insensitive to the value of r//s in the QGP 
phase but depends strongly on the value of rj/s in the hadronic phase. Our results show that, 
at least under the approximations made here, for the thermal photon spectrum the situation 
is quite the opposite: the spectrum is much more sensitive to the QGP shear viscosity than to 
the values in the hadronic phase. As shown in Figure HI the temperature dependence of the 
ratio rj/s in the hadronic phase (which corresponds to Model B) leads to a thermal photon 
spectrum that is almost indistinguishable from the spectrum computed from a temperature- 
independent value of t]/s (Model A), in sharp contrast to what happens with the spectrum 
obtained from an rj/s that varies in the QGP phase (which corresponds to Model C). One 
should be careful to note that small pr photons are mostly produced during the hadronic 
phase, whereas large px photons are mainly produced in the QGP phase. Therefore, it is 
perhaps more accurate to conclude that large px photons are more sensitive than small pr 
photons to variations in rj/s . This result implies that serious difficulties would arise when 
attempting to extract the temperature-dependence of r]/s corresponding to the hadronic 
phase by matching fluid dynamic simulations to measured photon spectra. 

Although it is tempting to conclude that fitting the spectrum of thermal photons to 
fluid dynamic simulations based on different parametrizations of temperature-dependent 
rj/s would allows us to directly extract this temperature dependence in the QGP phase 
from data, in reality much more work is needed to reach this stage, essentially because it 
is still necessary to determine more precisely the uncertainty coming from various sources 
in viscous fluid dynamic simulations. An example of such a source of uncertainty in the 
possible extraction of rj/s from the measured photon spectrum is presented later on when 
discussing results obtained with different EOSs. 

Returning to our results, Figure [3] shows the total thermal photon spectrum corresponding 
to the three models for rj/s. Also shown are data of direct photons as measured by the 



Phenix Collaboration 



241 ] . We note that we do not attempt to fit our results to data, since as 



mentioned in the Introduction we neglect prompt and jet-medium photons in calculating the 
total spectrum, considering only thermal photons. The rationale for showing experimental 
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results is merely to show that our viscous fluid dynamic model is not excluded by data, 
which would definitively be the case if the computed spectrum were larger that the observed 
one. It is seen from Figure E] that, as expected from the results shown in Figure [2J the 
total spectrum corresponding to Model C is considerably larger than those corresponding 
to Models A and B. Moreover, the difference between the results obtained from Model C 
and Models A-B increases with increasing values of pr and becomes significant for px > 1-5 
GeV. 

To put our results into perspective, it is useful to compare the impact that the 
temperature-dependence of r]/s has on thermal photon spectra to the effect of another in- 
put of fluid dynamic simulations, namely the EOS. This allows us to determine the relative 
importance of the temperature-dependence of rj/s on the spectra as compared to another 
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FIG. 2: (Color online) Contributions to the thermal photon spectrum corresponding to the QGP 
(upper panel) and the hadronic phases (lower panel), for Models A (temperature-independent), B 
(varying in the hadronic phase) and C (varying in the QGP phase) for rj/s. 
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important source of uncertainty in fluid dynamic simulations as is the precise shape of the 
EOS of nuclear matter created in heavy ion collisions. 

To this purpose, we present results of photon spectra obtained from fluid dynamic simula- 
tions that use two different EOSs as input. Figure J3] shows the spectrum of thermal photons 
obtained using the EOS of Laine and Schroder [43] (we will call it LQCD EOS in what 
follows) and the EOS shown in the lower panel of Figure HJ The latter EOS corresponds 



to the linear sigma model and was obtained from the linearized Boltzmann equation in 



23| 



we will call it EOS 1), and is used here for illustrative purposes. We note that EOS 1 



has been used in 



23| to study the influence of chiral fields on charged hadron observables. 



The upper panel of Figure H] shows the contribution of the QGP and hadronic phases to the 
total spectrum, which is shown in the lower panel. To keep things simple, we present the 
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FIG. 3: (Color online) Total thermal photon spectrum corresponding to Models A (temperature- 
independent), B (varying in the hadronic phase) and C (varying in the QGP phase) for rj/s. The 
data corresponds to direct photons as measured by the Phenix Collaboration 24 . 



11 



comparison for a temperature-independent rj/s — 0.11 (corresponding to Model A). 

It is seen that the main difference between the spectra obtained with the EOSs comes 
from the hadronic phase, with the photon spectrum being considerably smaller in the case of 
EOS 1. The difference between the spectra in both models is practically independent of pt, 
except at very low values where this difference becomes slightly larger. There are two reasons 
that make the spectrum obtained from the EOS 1 model smaller than the corresponding 
to the LQCD EOS. First, the hydrodynamic evolution is much faster in the EOS 1 model, 
essentially because the speed of sound is on average larger than the one corresponding to 
the LQCD EOS. This is illustrated by the fact that the freeze-out temperature is reached 
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FIG. 4: (Color online) Thermal photon spectrum obtained from simulations using the EOS of Laine 
and Schroder 43 (denoted by LQCD) and an EOS computed from a simple chiral-hydrodynamic 
model 23 (denoted by EOS 1). The upper panel shows the contributions corresponding to the QGP 
and hadronic phases, while the lower panel shows the total spectrum. The results correspond to a 
temperature-independent rj/s = 0.11 (Model A). 
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in ~ 6 fm/c in the EOS 1 model and in ~ 10 fm/c in the LQCD EOS, which represents 
a significant difference. The other reason is that the shear tensor U^ u is on average much 
smaller in the EOS 1 model, making the dissipative contribution to photon production 
smaller. Therefore, the spectrum of thermal photons is sensitive to the EOS used in the 
fluid dynamic simulations in two ways, namely through viscous corrections to the spectrum 
and through the rate of expansion and cooling of the fireball. 

Comparing the results for different parametrizations for the temperature dependence of 
rj/s with those obtained using different EOSs, it is seen that the effect of the temperature 
dependence of rj/s on thermal photon spectra becomes appreciable only at px > 1.5 GeV, 
whereas the effect of using different EOSs is significant in the whole range of values of pt in 
the case of the hadronic contribution and very small for the QGP contribution. Moreover, 
it is also seen that when the difference between the results obtained using Models A, B 
and C is appreciable, the influence on the spectra of the temperature dependence of rj/s 
and of the EOS is comparable. This implies that, if one would attempt to extract rj/s by 
matching the thermal photon spectrum computed from fluid dynamic simulations to data, 
the uncertainty associated to the temperature dependence of the ratio and to the EOS 
used in the simulations would be comparable. Since neither the EOS nor the temperature 
dependence of rj/s of nuclear matter created in heavy ion collisions is precisely known, it 
would prove quite difficult to extract the temperature dependence of rj/s from measured 
photon spectra. Our results highlight the necessity of carrying out a combined study of 
different observables (photonic and hadronic) in order to obtain more accurate constraints 
on the temperature dependence of the shear viscosity-to-entropy-ratio of matter created in 
heavy ion collisions. 

IV. CONCLUSIONS 

In this paper we have computed the spectrum of thermal photons including dissipative 
corrections in the production rates, using three parametrizations for the temperature depen- 
dence of rj/s that enter as input in fluid dynamic simulations that employ a realistic Lattice 
QCD equation of state. We have also compared the influence of the temperature depen- 
dence of rj/s on spectra to that of using different equations of state in the hydrodynamic 
simulations. 
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Our results suggest that some caution must be exercised when attempting to extract 
values of rj/s by fitting viscous fluid dynamic simulations to measured photon spectra. By 
using simple parametrizations for 77/s as a function of temperature, we have shown that 
the temperature dependence of this ratio has a significant influence on the thermal photon 
spectrum, which increases with increasing transverse momentum, becoming significant for 
Pt > 1.5 GeV. Probably the most important result that we obtain is that the spectrum of 
thermal photons depends significanlty on the values of rj/s in the QGP phase but is quite 
insensitive to the values in the hadronic phase, a feature that deserves further study. Work 
is in progress along this line. 
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